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Abstract 

We present a combined measurement of the Cabibbo-Kobayashi-Maskawa matrix element \V c b\ and 
of the parameters p 2 , Ri, and R2, which fully characterize the form factors of the B° — » D*~l + V£ 
decay in the framework of HQET, based on a sample of about 52,800 B° — > D*~£ + vi decays 
recorded by the BABAR detector. The kinematical information of the fully reconstructed decay 
is used to extract the following values for the parameters (where the first errors are statistical 
and the second systematic): p 2 = 1.156 ± 0.094 ± 0.028, i?i = 1.329 ± 0.131 ± 0.044, R 2 = 
0.859 ± 0.077 ± 0.022, ^(l)^ = (35.03 ± 0.39 ± 1.15) x 10~ 3 . By combining these measurements 
with the previous BABAR measurements of the form factors which employs a different technique on 
a partial sample of the data, we improve the statistical accuracy of the measurement, obtaining: 
p 2 = 1.179 ± 0.048 ± 0.028, R x = 1.417 ± 0.061 ± 0.044, i? 2 = 0.836 ±0.037 ±0.022, and ^(1)1^1 = 
(34.68 ±0.32 ±1.15) x 10~ 3 . Using the lattice calculations for the axial form factor ^(1), we extract 
\V cb \ = (37.74 ± 0.35 ± 1.25±};|1) x 10" 3 , where the third error is due to the uncertainty in ^(1). 
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1 INTRODUCTION 



The study of the B° — > D*~l + vp decay [1] is interesting in many respects. In the Standard Model 
of electroweak interactions the rate of this weak decay is proportional to the absolute square of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V c b, which measures the weak coupling 
of the b to the c quark. Therefore the determination of the branching ratio of this decay allows for 
a determination of \V c b\- 

The decay is also influenced by strong interactions, whose effect can be parameterized through 
two axial form factors A± and A2, and one vector form factor V, each of which depend on the 
momentum transfer q 2 of the B meson to the D* meson. The form of this dependence is not known 
a priori. Equivalently, instead of q 2 , the linearly related quantity w, which is the product of the 
four- velocities of the B° and D* + , defined in Eq. (1), can be used. 

In the heavy-quark effective field theory (HQET) framework [2, 3], these three form factors 
are related to each other through heavy quark symmetry (HQS), but HQET allows for three free 
parameters which must be determined by experiment. 

Several experiments have measured \V c b\ based on the study of the differential decay width 
dT/dw of B° -► D*-£+v e decays [4, 5, 6, 7, 8, 9, 10]. They determine only one of the parameters of 
the HQET from their data, while they rely on the CLEO measurement of the other two, as described 
in Ref. [11]. This is the largest systematic uncertainty common to all previous measurements of 
\Vcb\ based on this method. 

Improved measurements of all parameters are also important for describing the dominant back- 
ground from b — > c — > t cascades in both inclusive and exclusive b — > ulv decays for measurements 
of \V u b\. These motivations have been the basis of the measurement of all three parameters per- 
formed by BABAR and presented in Ref. [12]. While a new value of \V c j,\ has been computed using 
the new results as input to the previous measurement, the full correlation between \V c b\ and the 
three form factor parameters cannot be fully accounted for in that approach. 

In the present analysis, we perform a simultaneous measurement of both \V c b\ and of the form 
factor parameters from the measurement of the full four-dimensional differential decay rate (see 
Sec. 2). We extend the method of Ref. [10] by combining the measurement of different one- 
dimensional decay rates, fully accounting for the correlations between them, to extract the whole 
set of observables. 

A brief outline of the paper is as follows: In Sec. 2, we introduce the definition of all used 
observables, form factors, form factor combinations, and parameters. In Sec. 3 we briefly describe 
the relevant aspects of the detector and datasets we utilize in this measurement. In Sec. 4 we 
describe event reconstruction and selection and in Sec. 5 the analysis method. Sec. 6 describes fit 
results and the estimation of systematic uncertainties. In Sec. 7 we summarize the results of this 
analysis and combine form factor measurements with a previous BABAR measurement to arrive at 
a further improvement of the statistical errors, both for the form factor parameters and for the 
measurement of \V c b\. 

2 FORMALISM 

This section outlines the formalism and describes the parameterization used for the form factors. 
More details can be found in Refs. [2, 3]. The lowest order quark- level diagram for the decay 
B° — > D*~l + ve is shown in Fig. 1. 
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Figure 1: Quark-level diagram showing the weak interaction vertices in the decay B — ► D* £ + is£. 



2.1 Kinematic variables 

The decay B° — > D*~l + V£ is completely characterized by four variables, three angles and q 2 , the 
square of the momentum transfer from the B to the D* meson. 

The momentum transfer is linearly related to another Lorentz-invariant variable, called w, by 

Pb-Pd* m% + m 2 D , - q 2 

w = v B -v D * = = —2- , (1) 

rriBtriD* ZmBrriD* 

where and wifl. are the masses of the B and the D* mesons, pb and p&* are their four-momenta, 
and vb and vd* are their four- velocities. In the B rest frame the expression for w reduces to the 
Lorentz boost 7£>* = £x>*/Md». 

The ranges of w and q 2 are restricted by the kinematics of the decay, with q 2 ~ corresponding 

to 

m 2 D + mn* , n 

w max = / D « 1.504 2 
LrriB'mD* 

and w m i n = 1 corresponding to 

<La* = {™ B ~ m D *f « 10.69 ( GeV/c 2 ) 2 . (3) 




Figure 2: Kinematics of a B — > D* decay, mediated by an intermediate vector boson W; ir s 

is the soft pion from the D* decay. This diagram defines the three angles 6i, Qy and \- 

In this analysis we only reconstruct the decay D* + — > D°7r + , where — > i^7r, D° — > Kitttit or 
£)° — > Ktttt . The angular variables, shown in Fig. 2, are 

• 0£, the angle between the direction of the lepton in the virtual W rest frame and the direction 
of the virtual W in the B rest frame, 



• 9y, the angle between the direction of the D in the D* rest frame and the direction of the 
D* in the B rest frame, 

• X) the dihedral angle between the plane formed by the D* — D and the plane formed by the 
W — I system. 

2.2 Four-dimensional decay distribution 

The Lorentz structure of the B° — > D*~£ + i/£ decay amplitude can be expressed in terms of three 
helicity amplitudes (H + , and Hq), which correspond to the three polarization states of the 
D* (two transverse and one longitudinal). For light leptons these amplitudes are expressed [2, 3] 
in terms of the three form factors h Al , Ri and R2 ■ 

H i {w)=m B \ = ^h Al (w)Hi(w), (4) 

2v 1 — 2wr + r z 

where 



VI - 2wr + r 2 (l ± y^y^i 



= < 5 > 

1 — r 



2Jm B m D * m D * 
where it = — 1 and r — 



mB + m D * ms 
The full differential decay rate in terms of the three helicity amplitudes is [2, 3]: 

dT(B° D*-£+v e ) 6m B m 2 Dt rVw 2 - 1(1 - 2wr + r 2 ) 

dwd cos 6» £ dcos 6> v dx ~ 8(4vr) 4 ^ ^ 

x G||y c 5| 2 [(l -cos^) 2 sin 2 y iT 2 H 

+ (I + cos 6 e ) 2 sin 2 9 v H 2 (w) 

+ 4sin 2 6^cos 2 # V '#oM 

- 2sin 2 6 e sm 2 6 v cos2xH + (w)H„(w) 

— 4 sin Oe(l — cos 0^) sin#y cos Qy cos \ x H + (w)Hq(w) 

+ 4 sin 6^(1 + cos 6fy) sinfly cos 0y cos^ x H-(w)Hq(w)] , 

where all three of the Hi are functions of w. The four-dimensional distribution of w, cos 0£, cos 6>y, 
and x described by Eq. (7) is the physical observable from which we extract the form factors. The 
normalization of this distribution is directly related to \V c b\. 

2.3 Form- factors parameterizations 

The HQET does not predict the functional form of the form factors, and one needs a parameter- 
ization of these form factors. In the limit of infinite mass for the b and c quark, R\ and R2 are 
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simply both equal to 1. Corrections to this approximation are calculated in powers of (Aqcd/uiq) 
and as- 

Parameterizations as a power expansion in (w — 1) have been proposed in literature. The first 
linear parameterization is given by: 

h Al H = h M {l)(l- p 2 {w-l)). (8) 

A generic second order extension can be realized adding an extra parameter c as coefficient of 
the (w — l) 2 term: 

h Al (w) = h Al (l)(l- p 2 (w-l)+c(w-l) 2 ). (9) 

In both cases R\ and R2 are considered to be constants, neglecting their w dependence. 

Among various predictions relating the coefficients of higher order terms to that of the linear 
term, we decided to adopt the Caprini, Lellouch and Neubert [13] one (in the following referred to 
as CLN); the authors find the following functional forms for the three HQET form factors: 

h Al (w) = h Al (l)[l-p 2 hAi Z + (53p 2 hAi -15)z 2 (10) 

-(231^ -91)z 3 ], 
R^w) = R^l) - 0.12(it; - 1) + 0.05(w - l) 2 , 
R 2 (w) = i? 2 (l) + 0.11(tt;-l)-0.06(u;-l) 2 , 

where: 

V^TT-V2 

z = , 7= ■ (11) 

v^TT+v/2 

These functional forms are determined apart from three unknown parameters, p\ A , -Ri(l) and 
i?2(l), which must be extracted from the data. 

It is important to notice that h Al {\) = J-'(l), corresponding to the value of the Isgur-Wise 
function evaluated at zero recoil (using the notation commonly found in literature for dT/dw). 
A recent lattice calculation [14] (including a QED correction of 0.7%) gives: h Al (l) = = 
0Q1Q+ - 030 



3 THE BABAR DETECTOR 

The BABAR detector is described in detail in Ref. [15]. The momenta of charged particles are 
measured by a tracking system consisting of a five-layer silicon vertex tracker (SVT) and a 40- 
layer drift chamber (DCH), operating in a 1.5-T solenoidal magnetic field. Charged particles of 
different masses are distinguished by their energy loss in the tracking devices and by a ring-imaging 
Cherenkov detector. Electromagnetic showers from electrons and photons are measured in a CsI(Tl) 
calorimeter. Muons are identified in a set of resistive plate chambers inserted in the iron flux-return 
yoke of the magnet. 

The analysis is based on a data sample of 79 fb _1 recorded on the T(4S) resonance, and 9.6 ftT 1 
recorded 40 MeV below it, with the BAEL4R detector [15] at the PEP-II asymmetric-energy e + e~ 
collider. We use samples of GEANT Monte Carlo (MC) [16] simulated events that correspond to 
about three times the data sample size. 
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4 RECONSTRUCTION AND EVENT SELECTION 



The reconstruction of the events is the same as used in Ref. [10], and the event selection is an 
improved version of the one used in that paper. 

We select events that contain a D* + candidate and an oppositely charged electron or muon 
with momentum 1.2 < pe < 2.4 GeV/c. Unless explicitly stated otherwise, momenta are measured 
in the T(4S) rest frame, which does not coincide with the laboratory frame due to the boost of 
the PEP-II beams. In this momentum range, the electron (muon) efficiency is about 90% (60%) 
and the hadron misidentification rate is typically 0.2% (2.0%). We select D* + candidates in the 
momentum range 0.5 < po* < 2.5 GeV/c in the channel D* + — > D irf, with the D° decaying 
to K~it + , K~Tr + Tr~Tr + , or K~it + it® . The charged hadrons of the D candidate are fitted to a 
common vertex and the candidate is rejected if the fit probability is less than 0.1%. We require 
the invariant mass of the hadrons to be compatible with the D° mass within ±2.5 times the 
experimental resolution. This corresponds to a range of ±34MeV/c 2 for the D° — > K~ir + ir° decay 
and ±17MeV/c 2 for the other decays. For the decay D° — > K~tt + tt , we accept only candidates 
from portions of the Dalitz plot where the square of the decay amplitude, as determined by Ref. [17], 
is at least 10% of the maximum it attains anywhere in the plot. For the soft pion from the D* + 
decay, 7r+, the momentum in the laboratory frame must be less than 450 MeV/c, and the transverse 
momentum greater than 50 MeV/c. Finally, the lepton, the 7r+, and the D° are fitted to a common 
vertex using a beam-spot constraint, and the probability for this fit is required to exceed 1%. 

In semileptonic decays, the presence of an undetected neutrino complicates the separation of 
the signal from background. We compute a kinematic variable with considerable power to reject 
background by determining, for each 5-decay candidate, the cosine of the angle between the mo- 
mentum of the B° and of the D* + £~ pair, under the assumption that only a massless neutrino is 
missing: 

2E m E DH - M| - M 2 DH 

COS @B ,D*£ = o • 

This quantity constrains the direction of the B° to lie along a cone whose axis is the direction of 
the D* + £" pair, but with an undetermined azimuthal angle about the cone axis. The value of w 
varies with this azimuthal angle; we take the average of the minimum and maximum values as our 
estimator w for w. This results in a resolution of 0.04 on w. We divide the sample into 10 bins in 
w from 1.0 to 1.5, with the last bin extending to the kinematic limit of 1.504. 

The selected events are divided into six subsamples, corresponding to the two leptons and the 
three D° decay modes. In addition to signal events, each subsample contains background events 
from six different sources: 

• combinatoric (events from BB and continuum in which at least one of the hadrons assigned 
to the D* + does not originate from the D* + decay); 

• continuum (D* + £~ combinations from e + e~ — > cc); 

• fake leptons (combined with a true D* + ); 

• uncorrelated background (£ and D* + produced in the decay of two different B mesons); 

• events from charged B background, B + — > D*£ + v^X (via B + — ► D**^ i + f£ or non-resonant 
B + — > D*~t + h>£K + production), and neutral B background B° — > D*~£ + ^X, (via B° — > 
D**-g+ Vl or non-resonant B° — > D*'£ + U£7r° production); 
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140 150 160 

AM (MeV/c 2 ) 

Figure 3: Yields of on-resonance data (points) and the results of the fit (line) to the AM distri- 
bution, with contributions from continuum, fake-lepton, and combinatoric-D* + background events 
summed over all w bins. 



• and correlated background events due to the processes B° — > D* + ut , r — > £ X and 
B° _^ D*+X c , X c -► £~Y. 

We estimate the correlated background (which amounts to be less than 0.5% of the selected can- 
didates) from the Monte Carlo simulation based on measured branching fractions [19, 20], while 
we determine all others from data. Except for the combinatoric background, all other background 
sources exhibit a peak in the AM = M D *+ — M D o distribution, where M D *+ and M D o are the 
measured D* + and D° candidate masses. 

We determine the composition of the subsamples in each w bin in two steps. First we estimate 
the amount of combinatoric, continuum, and fake-lepton background by fitting the AM distri- 
butions in the range 0.139 < AM < 0.165 GeV/c 2 simultaneously to three sets of events: data 
recorded on resonance, data taken below the T(4S) (thus containing only continuum background), 
and data in which tracks that fail very loose lepton-selection criteria are taken as surrogates for 
fake leptons. The distributions are fitted with the sum of three Gaussian functions with a common 
mean and different widths to describe D* + — > D°irf decays and with empirical functions 

F comb (AM) = 1 

based on simulation for the combinatoric background. The four parameters of the Gaussian func- 
tions are common, while the fraction of peaking background events and the parameters describing 
the combinatoric background differ for the signal, off-peak, and fake-lepton samples. 

Since the AM resolution depends on whether or not the 7r+ track is reconstructed only in the 
SVT or in the SVT and DCH, the fits are performed separately for these two classes of events. 
We rescale the number of continuum and fake-lepton events in the mass range 0.143 < AM < 
0.148 GeV/c 2 , based on the relative on- and off-resonance luminosity and measured hadron misiden- 
tification probabilities. The fit is performed in multiple steps: the combinatoric background is first 
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Figure 4: Yields of on-resonance data (points) and the results of the fit (histograms) to the 
cos 9 bo d*£ distribution, summed over all w bins. 



fixed to the simulation prediction in order to optimize the fit of the signal in the tail region around 
the peak. Then the peak shape is fixed to the fitted one, and the combinatoric background is left 
free to float, to allow for possible small deviations from the simulation. 

In the subsequent analysis we fix the fraction of combinatoric, fake-lepton, and continuum 
events in each w bin to the values obtained following the procedure described above. Figure 3 
shows the AM fit results for the on-resonance data. 

In a second step, we fit the cos9 B o D *£ distributions in the range —10 < cos9 B o D * e < 5 and 
determine the signal contribution and the normalization of the uncorrelated and B — ► D*£ + V£X 
background events. Neglecting resolution effects, signal events meet the constraint | cos 9 B o £>*^| < 1, 
while the distribution of B — > D*l + vnX events extends below —1, and the distribution of the 
uncorrelated background events is spread over the entire range considered. 

We perform the fit separately for each w bin, with the individual shapes for the signal and for 
each of the six background sources taken from MC simulation, specific for each of the six subsamples. 
Signal events are generated with the CLN form factor parameterization, tuned to the results from 
CLEO [11]. Radiative decays (B° -► D* + £~ua) are modeled by PHOTOS [22] and are treated 
as signal. B — ► D**£v decays involving orbitally excited charm mesons are generated according 
to the ISGW2 model [24], and decays with non-resonant charm states are generated following the 
prescription in Ref. [25]. To reduce the sensitivity to statistical fluctuations we require that the 
ratio of B — > D*l + u^X and of uncorrelated background events to the signal be the same for all 
three D° decay modes and for the electron and muon samples. Fit results are shown in Fig. 4. In 
total, there are 68,840 events in the range | cos 9 B o D *(\ < 1.2. The average fraction of signal events 
is (76.7 ± 0.3)%, where the error is only statistical. 
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5 THE EXTRACTION OF | V c6 | AND THE FORM FACTOR PA- 
RAMETERS 



The approach followed in this analysis to simultaneously measure \V c b\ and the three form factor 
parameters is to extend the one-dimensional fit procedure of Ref. [10] to several one-dimensional 
distributions, while fully accounting for the mutual correlations, using the same sample of selected 
events. This permits disentangling the effects coming from each of the parameters, and to measure 
all of them. 

In principle, the multi-dimensional phase space could be divided into several independent bins, 
and the fraction of signal events in each of them could be determined as described in the previous 
section. In practice, the need for a detailed bin segmentation would imply serious statistical lim- 
itations in most of the bins. For this reason, we instead consider projections of the data into the 
most significant variables ( w, cos0£, cos#y). We divide each projection into ten bins, and deter- 
mine the amount of signal events in each bin applying the procedure described above. We then 
fit the resulting projections to extract the form factors and \V c b\. We account for the correlations 
between bins of different projections by noting that the covariance between the content of two bins 
in different one-dimensional distributions is given by the common number of events in these bins, 
while it is zero for bins in the same distribution. 

Any kinematic observable sensitive to the form factor parameters can be used, the sensitivity 
to \Vcb\ being already present in the w distribution alone. We have studied the distributions of 
the kinematic observables describing the decay, as discussed in Sec. 2. We have verified that the 
distribution of x is practically insensitive to the form factor parameter changes, therefore we have 
chosen w, cos0£ and cos 0y as set of observables to be included in the fit. Figure 5 shows the 
sensitivity of each of the 4 kinematical variables describing the B° — > D*~l + V£ decay to the form 
factor parameters. 

The value of \V c b\ and the form factor parameters are extracted from a least squares fit to the 
set of measurements given by the measured contents of bins of w, cos0£ and cos 6>y, according to 
the selection described in Sec. 4 and the cut | cos##o jj*g\ < 1.2. The w and cos0y distributions are 
divided into ten equal size bins. The cos0£ distribution is divided in the analysis in 10 non-equal 
size bins, whose boundaries are [—1.0, —0.4, —0.2, 0.0, 0.2, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0], chosen according 
to the strong variations in the shape, and therefore to the available statistics. 

5.1 The fit method 

The x 2 function to be minimized is an extension of the one used in Ref. [10]. Defining n as the total 
number of observable bins, Nf ata the total number of observed events in bin i, N^ k9 the number of 
estimated background events for that bin, and N^ 1 the number of simulated signal events falling 
in that bin, the x 2 can be written as: 




N MC 



(13) 




N MC 



15 





Figure 5: Sensitivity of the w, cos Op, cosOy and \ variables to the variations of the three form 
factor parameters p 2 , R\ and i?2- The variations considered are 1 a for p 2 , Ri, and R2 using the 
measured value from [10] for p 2 and from Ref. [11] for the others. 
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where W' is the weight assigned to the k-th simulated signal event falling into bin i in order to 
evaluate the expected signal yield as a function of the searched parameters, and Cij is the covariance 
matrix element for the bin pair ij. 

Each weight W k is the product of four weights, W k = W £ W^ ,fc wf ,fc wf /,fc . The factor W £ 
accounts for the normalization of data and simulation samples. It depends on the total number of 
BB events, N BB = (85.9 ± 0.9) x 10 6 , on the fraction of B°B° events, fa = 0.489 ± 0.012 [19], 
on the branching fraction B{D* + -► D°ir + ) = 0.677 ± 0.005 [19], and on the B° lifetime t b o = 
1.532 ± 0.009 ps [19]. The factor W^' fc accounts for differences in reconstruction and particle- 
identification efficiencies predicted by the Monte Carlo simulation and measured with data, as a 

function of particle momentum and polar angle. 

S k 

The factor W i ' accounts for potential small residual differences in efficiencies for the six data 
subsamples (3 D° decay modes and 2 leptons) and allows for adjustment of the D° branching 
fractions, properly dealing with the correlated systematic uncertainty. It is the product of several 
scale factors that are floating parameters in the fit, each constrained to an expected value with 
a corresponding experimental error. To account for the uncertainty in the multiplicity-dependent 
tracking efficiency, we introduce a factor Wj^ fc = 1 + Nt r k(^ — Strk), where Nt r k is the number of 
charged tracks in the D*£ candidates in each sample and 5t r -k is constrained to one within the esti- 
mated uncertainty in the single-track efficiency of 0.8%. Similarly, multiplicative correction factors 
are introduced to adjust lepton (SpiDe, fipiDfi) efficiencies, kaon efficiencies in the different selec- 
tions used for different D° decay channels (Spm Ktight, ^pid Knot apion), and ir° (5 n o) efficiencies, 
within the total estimated uncertainties of these efficiencies, and D° branching fractions ($brd° Km 
$brd° Knmn $br D° Kwn° ) ■ Correlations between the branching fractions are taken into account in 
the constraint through the covariance matrix Cbr (given by the uncertainties and the correlation 
matrix obtained from Ref. [19, 20]). Therefore wf' fc can be written as: 

Wf k = 8p IDe 5 PID ^ P k IDK {l + N l t ; k k (l - 5 trk ))6$5 B k RD0 (14) 

where the kaon and 7r° efficiency corrections, the decay multiplicity N trk and the D° branching ratio 
correction depend on the particular event k of the bin % considered (because of the different decay 
multiplicity, branching ratio, kaon selector and presence of 7r°). The addition of the constraints 
effectively adds extra terms to the above x 2 expression for each of the corrective factors S. The 
complete \ 2 function used in the fit therefore is: 

x 2 = N ? ata - N i kg - E w x c ^ ( N t ata - N f g - 11 w i ) ( 15 ) 

i=l j=l ^ k=l ' ^ k=l ' 

(l-5 PIDe ) 2 (l-5 PIDll ) 2 

2 2 
a PIDe (J PID i i 

(1 - & P ID K tight) 2 (1 - 5 pip 

K not apion ) 



+ 



2 2 
® PID K tight a PIDK not apion 

(l-5 trk ) 2 , (1-M 2 



+1 — + 



a lk a lo 



+ ^2 S BRD°m>< C BRD°mn S BRD°n 

m=K-K,K-K7nr,K-KiT n=Kir,Kmrir,KTrir 

The addition of these extra terms allows us to fit all subsamples simultaneously fully taking into ac- 
count the correlated systematic uncertainties, and effectively propagates these uncertainties through 
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the weights to the uncertainty on the free parameters of the fit. 

The fourth factor w{^' k accounts for the dependence of the signal yield on the parameters to be 
fitted. The \V c b\ dependence is trivially given by the ratio | V^ 2 /! V^f c \ 2 , where the denominator 
is the actual value used in the simulation, derived from the branching ratio used for the decay 
B° — ► D*~£ + U£. In an analogous way the dependence on the form factor parameters (p 2 , Ri, and 
R2) is given by the ratio of the differential decay rate of Eq. (7) evaluated at the variable values 
(in the chosen fitting model) and at the simulation values (where the linear model was used). 

One important point is that not all bins of all used observables can be used at the same time: in 
fact, once all bins belonging to one observable are used, the overall normalization is fixed, therefore 
for the other observables one of the bins has to be dropped, since its content and covariance matrix 
part can be written as linear combination of all others. The choice of the bins to be dropped is 
arbitrary, and it has been verified that by dropping different bins exactly the same numerical results 
were obtained (within the numerical precision of the fit). 

Summarizing, three observables are used in the fit, w, cos0£, and cos 9y, whose ranges are 
divided into ten bins. One bin is excluded from the fit for the last two observables, for a total of 
28 bins used in the x 2 - 

This fit procedure has been tested on a signal toy Monte Carlo, and no evidence for systematic 
biases has been observed. 

5.2 The observables' covariance matrix 

A direct consequence of this method is that the covariance matrix for the measurements corre- 
sponding to the bins used is not diagonal, since the bins are not all statistically independent. The 
total covariance matrix is the sum of 3 separate matrices: one for the measured data yields, one 
for the estimated signal yield, and one for the evaluated background. The diagonal elements of the 
matrices are given by the uncertainty on the bin content itself: the covariance of bins belonging to 
the same observable is zero, the covariance of bins from different observables is the variance of the 
number of events that is common between the bins. 

The observed data matrix is built assuming the bin contents and their intersections obeying a 
Poisson distribution, therefore the variance of a bin content or of the content of the intersection 
of two bins is the content itself. The estimated signal matrix is built in an analogous way, where 
the variance of a number of weighted events n is approximated by the sum of the squared weights 

£<=!,„ w l 

The calculation of the background covariance matrix is less straightforward. The diagonal part 
is simply the estimated variance of the measured background, according to the whole procedure 
described in Sec. 4. The background extraction procedure does not directly determine the common 
number of events between two bins, because this procedure is based on a complex sequence of 
shapes and yield fits. 

The solution adopted is to use the number of common background events between two bins 
as predicted by the simulation, after it has been corrected for tracking and PID efficiencies, and 
rescaled in such a way as to get for each bin of one of the studied kinematic observables (one among 
w, cos 0£, cos 0y, x) the background amount estimated from the data-based background evaluation. 
In this estimation procedure also the x projection is used, in order to provide an additional check 
on the background evaluated from the three distributions used in the fit. 

Since the total integral of the background measured is not exactly the same when evaluated 
for different kinematic observables, it has to be fixed to a given value. The average normalization 
determined from the four kinematic observables describing the event is used. 
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The spread of the background normalization values is found to be about two times larger than 
the estimated uncertainty from the background procedure error propagation. This is due to the fact 
that the uncertainties on the signal and background shapes are not accounted for in the background 
estimate. Taking into account the structure of the distribution of the four integrals, split into two 
subgroups of compatible values, the global background uncertainty is estimated to be the average 
of the maximum and minimum difference between these two groups (i.e., w and cos6*£ on one side, 
cos 6y and x on the other) . The background uncertainties are suitably rescaled in order to give the 
uncertainty on the integral equal to the one estimated as explained. 

6 THE FIT RESULTS AND SYSTEMATIC UNCERTAINTIES 

6.1 Results of the fit 

The analysis result is obtained using the background covariance matrix evaluated from the measure- 
ment in the variable w. Table 1 shows the best fit result. The uncertainties have been determined 
by the MINDS [21] algorithm, whose output has been symmetrized. The 5 coefficients describing 
systematic uncertainties not common to all events are found to be consistent with 1 within the 
estimated uncertainties. The contribution to the overall uncertainty due to statistical fluctuations 
of the data and simulated samples has been obtained by recomputing the covariance matrix of the 
fit after fixing all the 5 parameters at their best fit values. 

Figure 6 shows the distributions of the kinematical observables describing the B° — > D*~£ + i/£ 
decay with the simulation prediction computed using the best fit values for \V c b\ and for the form 
factors. 

The B° — ► D*~l + i/£ decay branching ratio is shown in the last row of the Table 1. The 
uncertainties are obtained through the propagation of the ones on p ,R±, R2 and ^"(1)1^1, taking 
into account their correlations. 



Parameter 


Best value 


<?tot 


O'stat 


&syst 


P 2 


1.156 


0.095 


0.094 


0.004 


Ri 


1.329 


0.131 


0.131 


0.007 


R2 


0.859 


0.077 


0.077 


0.002 


F(l)\V cb \ x 10 3 


35.03 


0.95 


0.39 


0.86 


BR(B° -► D*-£+vi) % 


4.83 


0.24 


0.05 


0.24 



Table 1: Fit result for the standard analysis. a to t is the total uncertainty, and the last two columns 
show its splitting in statistical and systematic contributions (from 5 parameters). In the last row 
the results for the branching ratio of B° — ► D*~£ + u^ obtained from the fit results are shown. 

Table 2 shows the correlation matrix for the standard analysis result. 

The goodness of fit can be evaluated directly from the value of the fit function in the minimum, 
since it is based on a least squares method, obtaining a x 2 /d.o.f. = 23.99/24 corresponding to a 
X 2 probability of 46.2%. In this evaluation the number of degrees of freedom is considered equal 
to the number of used bins, ignoring the bin correlations, minus the number of fitted parameters. 
It is interesting also to check the goodness of the fit on the separate distributions entering the fit 
(i.e., w, cosdi and cos 6>y), considering for each of them all ten bins in which their ranges have been 
subdivided. The corresponding results are given in Table 3. 
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Figure 6: Distributions of w, cos9i, cosOy, x after the fit. The error bars represent the measured 
values, the colored histograms represent the background estimate and the fitted signal MC. The 
bottom plots show the ratio between the measurement and the sum of all background and fitted 
signal MC events. 
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o 2 
P 




ll 2 


T(lMVu\ x 10 3 


P 2 


+1.000 


+0.864 


-0.922 


-0.028 


Ri 


+0.864 


+1.000 


-0.928 


-0.226 


R2 


-0.922 


-0.928 


+1.000 


+0.172 


F{l)\V cb \ x 10 3 


-0.028 


-0.226 


+0.172 


+ 1.000 



Table 2: Total correlation matrix for the standard analysis fit. 



Variable 


xVd.o.f. 


X 2 probability 


w 


3.62/10 


96.3% 


cos 9 1 


12.93/10 


22.8% 


COS 9y 


6.04/10 


81.2% 



Table 3: Goodness of fit for the distributions entering the fit for the standard analysis. All ten bins 
are considered. 

In Fig. 7 the contour plot at the 39% C.L. (A X 2 = 1) in the p 2 -Ri, p 2 -R 2 , R1-R2, p 2 -F{l)\V ch \ x 
10 3 , Vcb\ x 10 3 , ^2-^(1)1^1 x 10 3 planes are shown, illustrating the correlations among 

the variables. 



6.2 Systematic uncertainties 

The summary of systematic uncertainties on the measured parameters is presented in Table 4. In 
the last column, the uncertainty breakdown for the B° — > D*~l + vn branching ratio is also shown. 



CLN 


P 2 


Ri 


R2 


F(l)\V cb \ x 10 3 


BR(B° -► D*~e+v e ) 


stat. error 


0.094 


0.131 


0.077 


0.39 


0.05 


PID, trk, BR D u 


0.004 


0.007 


0.002 


0.86 


0.24 


7T Soft 


0.013 


0.005 


0.001 


0.46 


0.18 


D*l vertex 


0.014 


0.010 


0.008 


0.06 


0.06 


B mom. 


0.013 


0.040 


0.017 


0.29 


0.14 


rad. corr. 


0.005 


0.004 


0.000 


0.19 


0.07 


D** compos. 


0.011 


0.008 


0.009 


0.10 


0.07 


background cov. matrix 


0.006 


0.004 


0.002 


0.04 


0.04 


partial sum 


0.027 


0.043 


0.021 


1.04 


0.35 


B u lifetime 








0.10 


0.03 


BB number 








0.19 


0.05 


BR(D* -» D°tt) 








0.13 


0.04 


/++//00 


0.008 


0.010 


0.006 


0.41 


0.15 


tot. syst. error 


0.028 


0.044 


0.022 


1.15 


0.39 



Table 4: Breakdown of systematics uncertainties. 
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6.2.1 Statistical, PID and tracking efficiency and D branching fraction uncertainties 

The uncertainty on the parameters given by the fit is not of purely statistical origin, since the sys- 
tematic uncertainty sources which are not common to all events (PID efficiency, tracking efficiency, 
and D° branching fraction uncertainties) are taken directly into account in the fit through the 5 
coefficients appearing in the weights Wf' fc . The splitting between statistical and systematic part 
of the fit result uncertainty is obtained by computing the purely statistical covariance matrix (by 
fixing all 5 parameters at their best value) and subtracting it from the global one. 

6.2.2 Soft pion efficiency 

The uncertainty coming from the soft pion efficiency correction is determined by propagating the 
uncertainties of the parameters of the efficiency correction function. In addition, an overall efficiency 
uncertainty due to track reconstruction and selection is considered (a scale factor affecting only 
.F(l)|Vd>|), hs magnitude being 1.3%. The two uncertainties are added in quadrature. 

6.2.3 D*l vertex reconstruction efficiency 

The uncertainty coming from the vertex reconstruction method used to select the D*l system 
has been evaluated using two alternative choices to the standard method: the first corresponds to 
removing the beam spot constraint from the vertex fit, while the second corresponds to removing the 
lepton information from the vertex fit. The biggest variation is taken as the systematic uncertainty. 

6.2.4 B momentum 

A correction to the B momentum in the simulation has to be applied to obtain agreement with 
data. Effectively this implies rescaling the cos0b°,d*£ value determined for the simulation by a 0.97 
factor. Half of the relative effect of the correction is considered as systematic uncertainty on it. 

6.2.5 Radiative corrections 

Radiative corrections to the B° — ► D*~£ + vg decay are computed in the simulation by PHOTOS [22], 
which describes the final state photon radiation (FSR) up to 0(a 2 ). In the event reconstruction no 
attempt to recover photons emitted in the decay is performed, therefore the simulated prediction 
of the reconstructed event is sensitive to the detail of the radiative corrections. This is particularly 
important in the cos6 B o D *£ fit for electrons, where the FSR is responsible for the long tail at low 
values. 

No detailed calculation of the full 0(a) radiative corrections to the B° — > D*~£ + U£ decay is 
available in literature. Recently, a new 0(a) calculation of radiative corrections in kaon decays 
has become available [23] explicitly taking into account radiated photons, and a new detailed 
comparison with PHOTOS has been presented. From this comparison it is evident that the radiated 
photon's energy spectrum is quite well reproduced by PHOTOS, the main difference being in the 
angular distribution of the photons with respect to the lepton. 

In order to assess the systematic uncertainty due to the imperfect treatment of the radiative 
corrections, we have used the comparison presented in Ref. [23] to reweight our simulated events 
in order to reproduce the photon angular spectrum given there, for photons above 10 MeV in the 
B center of mass. The effect of the reweighting has been used as an estimate of the systematic 
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uncertainty. In this approach the possible effect of the B° decay form factors on the FSR description 
is completely ignored. 

6.2.6 B — ► D*t r viX background description 

The shapes of the different components of the B — ► D*l + uiX background are taken from the 
simulation. In order to account for the uncertainty on their knowledge, the fit has been performed 
using only one of them at a time. The study has been performed for the D° — > Kir subsample only, 
and assumed to be valid for all subsamples. The uncertainty on the result is given by half of the 
biggest variation observed. 

6.2.7 Fit covariance matrix 

The effect of using the background shape evaluated in different projections to compute the back- 
ground covariance matrix is used as an evaluation of the corresponding systematic uncertainty (the 
background overall normalization is fixed to the average of the one computed in the projections). 
The maximal variation with respect to the basic result when changing the projection used is taken 
as an estimate. 

6.2.8 Global normalization factors 

In the second part of Table 4 the effect of the uncertainty on quantities which are normalization 
factors is shown. The effect of the uncertainties in the D° branching ratios has already been dis- 
cussed, and it affects all parameters because it is changing the relative amount of signal events 
coming from different D° decays. On the other hand, the B° lifetime, the D* — > D°7r branching ra- 
tio, and the BB measured number are global scale factors affecting only ^"(l)|Vrf,|. The uncertainty 
on the ratio /++//oo = B(T{AS) -► B+ B~) /B(T{AS) -► B°B°) affects both the absolute number 
of measured B° — > D*~£ + U£ decays (and therefore JF(1)| V^l) and the relative ratio of background 
events from B° and B^ in the simulation. This second aspect influences the cos^go d*i shapes 
and therefore the background determination. For this reason this uncertainty source cannot be 
considered a pure scale uncertainty affecting only .F(l)| V c b\, but it touches also the form factors 
parameters. The analysis has been repeated by changing the /++//oo value by one sigma and the 
variation has been used as systematic uncertainty on the results from this source. The B° lifetime 
used is t b o = 1.532 ±0.009 [19]. The D* -► D°ir branching ratio used is 67.7 ±0.5%. The _/++//oo 
value used is 1.055 ±0.055 from Ref. [18]. The systematic uncertainty on the number of BB (N BB ) 
is 1.1%. The uncertainties on these input quantities are propagated to the final fit results. 

7 CONCLUSIONS 

7.1 Summary of results of this analysis 

A sample of about 52800 fully reconstructed B° — > D*~£ + U£ decays collected by the BABAR detector 
has been used to measure both ^(l^V^I and the form factor parameters in the CLN parameter- 
ization [13] p 2 , R\, and R2 in a global fit. Both identified electrons and muons have been used, 
the D*~ candidates are reconstructed from the D*~ — ► D° ir~ decay and the D° candidate is 
reconstructed in three different decay modes K~ir + , K~ir + ir + Tr~ and K~tt + it . 
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The obtained results are (where the first errors are statistical and the second systematic) : 



F{l)\V ch \ = (35.03 ± 0.39 ± 1.15) x 10~ 3 

p 2 = 1.156 ±0.094 ±0.028 

i?i = 1.329 ±0.131 ±0.044 

R 2 = 0.859 ±0.077 ±0.022 

The correlations between the fitted parameters are: 

p(p 2 ,Ri) = ±86% 

p{p\R 2 ) = -92% 

p(p 2 ,F(l)\V cb \) = -3% 

p(Ri,R 2 ) = -92% 

= -23% 

p(R 2 ,F(l)\V cb \) = ±17% 

Using a recent lattice calculation [14] (h Al (l) = JF(1) = 0.919t°;o35) results in the following 
value for \V cb \: 

\V cb \ = (38.12 ± 0.42 ± l.25+l 2 4 j) x 10~ 3 , 

where the third error is due to the uncertainty in J-{1). 

The corresponding branching fraction for the decay B° — > D*~£ + V£ is found to be: 

B(B° -» D*~£ + v e ) = (4.84 ± 0.05 ± 0.39)%. 



7.2 Combination of results with alternate form factor parameters measurement 

The BABAR collaboration recently published a measurement [12] of the same form factor parameters 
in B° — > D*~£ + V£ decays based on an unbinned maximum likelihood fit to the four-dimensional 
decay distribution using a subset of the data analyzed in this work. This earlier analysis has 
a higher sensitivity to the form factor parameters and resulted in p 2 = 1.145 ± 0.066 ± 0.035, 
Ri = 1.396 ± 0.070 ± 0.027, and R 2 = 0.885 ± 0.046 ± 0.013. 

We can combine the two form factor parameters measurements taking into account the corre- 
lation between them and obtain 



jc-(l)|y c6 | = (34.68 ± 0.32 ± 1.15) x 10~ 3 

p 2 = 1.179 ±0.048 ±0.028 

Ri = 1.417 ±0.061 ±0.044 

R 2 = 0.836 ±0.037 ±0.022. 

The statistical errors are significantly improved compared to the analysis presented in this paper. 
The two analyses have largely the same sources of systematic uncertainties, and thus we retain the 
systematic measurement errors established in this paper. The combined statistical errors are still 
larger than the systematic errors, but not by a large factor. 
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The correlation coefficients between the fitted parameters are: 



p{p\Ri) 
p(p 2 ,R2) 



+70% 
-83% 
+27% 
-84% 
-39% 
+22% 



p(p 2 ,F(l)\V cb \) 

p(Ri,R2) 
p(R 2 ,F(l)\V cb \) 



Using the lattice calculation for JF(1), we obtain an improved value for \V c b\ [14]: 



where the third error reflects the current uncertainty on J-(l). This combined measurement of 
\V c b\ represents a significant improvement over previous measurements. It supersedes all previous 
exclusive BABAR measurements based on partial data sets or less detailed analyses. 
The corresponding branching fraction of the decay B° — > D*~l + V£ is found to be: 
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